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Abstract 

Thermal inflation is an attractive idea to dilute cosmic density of unwanted 
particles such as moduli fields which cause cosmological difficulties. However, it 
also dilutes preexisting baryon asymmetry and some viable baryogenesis is necessary 
for a cosmologically consistent scenario. We investigate whether the Affleck-Dine 
mechanism can produce baryon asymmetry enough to survive after the dilution 
in gauge-mediated SUSY breaking models. Flat directions except for LH U flat 
direction cannot provide such huge baryon number because of Q-ball formation. 
We show that although the LH U flat direction is special in terms of having /r-terrn 
which prevents Q-ball formation, it cannot explain the observed baryon asymmetry 
either. 


1 Introduction 


Supersymmetry (SUSY) [1] is one of the most attractive candidates for extensions of the 
Standard Model, ft can not only naturally solve the hierarchy problem but also achieve the 
unification of the gauge couplings, which implies the existence of Grand Unified Theories 
(GUTs). Because no supersymmetric partners have been discovered so far, however, SUSY 
must be spontaneously broken, which generally causes a SUSY flavor problem. Gauge- 
mediated SUSY breaking models [2, 3, 4] give a solution to the flavor problem since the 
mediation by the gauge interactions leads to the degeneracy of sfermion masses. We focus 
on the gauge-mediated SUSY breaking models in this paper. 

In terms of the cosmology of SUSY models, some long-lived particles may cause serious 
problems. Local SUSY predicts the existence of gravitino, which is the superpartner of 
the graviton. Although initial abundance of gravitinos is diluted during the primordial 
inflation, they are produced again through scattering process and decay from sparticles 
during reheating. In gauge-mediation models, the gravitino mass m 3 / 2 is much lighter 
than the electroweak scale {m 3/2 0(100) GeV). It interacts with other particles only 
through gravitational interaction and hence has a long lifetime. I 11 particular, if the 
lifetime of the gravitino is larger than the present age of the universe, which is actually 
the case for m 3 / 2 < 0(100) MeV, its abundance may give too much contribution to the 
present cosmic density of the universe [5, 6]. In order to avoid this problem called “the 
gravitino problem”, the reheating temperature should be low enough since the produced 
gravitino abundance is proportional to the reheating temperature. 

In the framework of superstring theories, the situation gets worse. There exists a 
lot of singlet scalar fields corresponding to flat directions in superstring theories. We 
call them “moduli fields”. The moduli fields generally acquire mass comparable to the 
gravitino mass through some non-perturbative effects of the SUSY breaking [7]. They also 
interact with other particles only through gravitation (i.e., Planck suppressed operators) 
and hence have long lifetimes. Let us focus on one of the moduli fields and consider its 
dynamics in the early universe. During and after the inflation, vacuum energy which 
causes the expansion of the universe largely breaks SUSY and lifts the moduli potential. 
Then, the moduli field acquires mass of the order of the Hubble parameter and sits down 
at the minimum of the potential which is in general different from the one at present. 
When the Hubble parameter becomes of the order of the moduli mass, the moduli field 
starts to oscillate around the true minimum. Its initial amplitude is expected to be of 
the order of the Planck scale because the gravitational scale is the only scale appearing in 
supergravity actions. A large number of moduli particles are then produced as coherent 
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oscillations. The moduli field oscillates until present and their energy density leads to 
overclosure in the case of m v < 0(100) MeV, where m v is moduli mass. If the moduli 
field has mass of m v ~ 0.1 MeV-1 GeV, it generally decays into X-rays during the present 
epoch and it may give too much contribution to the X-ray background spectrum, which 
more severely constrains the abundance of the moduli [8, 9]. This cosmological difficulty 
is sometimes called “the cosmological moduli problem” [7, 10, 11]. An important point 
is that the moduli problem cannot be solved even if the reheating temperature is as low 
as 0(10) MeV, which is restricted from below to realize the Big Bang nucleosynthesis 
(BBN). The moduli problem is therefore more serious than the gravitino problem. 

In order to solve the moduli problem, late-time entropy production is needed to dilute 
the moduli abundance. A probable candidate for a huge entropy dilution mechanism is 
“thermal inflation” [12, 13] which is mini-inflation caused by a scalar field called “a flaton”. 
Introducing the flaton is also motivated in light of a solution to the fi problem [14, 15]: 
why is the /i -term electroweak scale and so small compared to the GUT or Planck scale? 
After the thermal inflation, flaton decay can produce a large amount of entropy enough to 
solve the moduli problem. Indeed, it has been studied that some specific models succeed 
in solving the moduli problem [16, 17]. 

However, one is faced with another problem incidental to the entropy dilution. The 
entropy production dilutes not only the moduli abundance but also baryon number which 
may be produced beforehand by some mechanisms. In particular, the thermal leptogen- 
esis [18, 19] cannot produce such huge baryon number that it can explain the observed 
baryon asymmetry after the dilution. Hence, we need a more effective baryogenesis or 
a mechanism to generate baryon number after the entropy production. One of the most 
probable candidates for the former mechanisms is the Affleck-Dine mechanism [20, 21], 
which can produce huge baryon asymmetry via dynamics of a flat direction in the visible 
sector. If a flat direction carrying baryon charge, called “the Affleck-Dine field”, acquires 
a nonzero vacuum expectation value (VEV) in the early universe, it starts to oscillate after 
the primordial inflation and rotates in the complex plane, which represents the generation 
of baryon number. 

The baryon number production by the Affleck-Dine mechanism is closely related to 
Q-ball [22] formation. Q-ball is a non-topological soliton that may be formed during the 
oscillation of the Affleck-Dine field [23, 24, 25]. It is well known that the Affleck-Dine field 
fragments into Q-balls if its potential becomes flatter than a quadratic term for a larger 
VEV. Since Q-balls absorb (almost) all produced baryon charge [25, 26, 27], we need to 
calculate baryon number released from these Q-balls to estimate the baryon asymmetry 
in the quark sector. 
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In this paper, we study if the Affleck-Dine baryogenesis is compatible with the thermal 
inflation in gauge-mediated SUSY breaking models. The potentials of the Affleck-Dine 
fields except for the LH U flat direction become flat at large VEVs in gauge-mediation 
models [6], which implies that Q-balls are formed after the Affleck-Dine baryogenesis. 
Unfortunately, it has been revealed that Q-ball formation rendered the Affleck-Dine mech¬ 
anism very unlikely to work in the presence of the moduli problem [28]. In this paper, 
we refine the estimated baryon number in [28] using a new lower bound of the SUSY 
breaking scale, and show that this scenario is more unlikely to work. On the other hand, 
the LH U flat direction is special in terms of having the /i-term which prevents it from 
forming into Q-balls. Although one can easily find that in the simplest case the resulting 
baryon asymmetry is not much enough to explain the observed abundance, we provide 
another scenario that provides larger amount of baryon asymmetry. Motivated to solve 
the /i problem [14, 15], we assume that the flaton couples with Higgs supermultiplets and 
its VEV generates the /i-term at low energy. In this case, the LH U flat direction starts 
to oscillate by a potential without the /i-term, which implies that it forms into Q-balls. 
Then, after the thermal inflation ends, the /i-term is generated and the LH U flat direction 
obtains the /i-term, which results in Q-ball deformation. We study if the released lepton 
number can explain the present baryon asymmetry in this scenario. 

The organization of this paper is as follows. In Sec. 2, we briefly explain the moduli 
problem and the thermal inflation. First, we consider a model of the thermal inflation 
introducing a flaton supermultiplet with a discrete symmetry. In this model, however, 
flaton decay into long-lived A-axions could lead to a cosmological difficulty in gauge- 
mediated SUSY breaking models. In order to prohibit the decay channel into R- axions, 
we introduce a symmetry breaking term as proposed in [16, 17]. In Sec. 3, we explain the 
Affleck-Dine mechanism as a candidate for an efficient production mechanism of baryon 
asymmetry. We also review Q-ball and its properties. In Sec. 4, we show that the Affleck- 
Dine baryogenesis is incompatible with the moduli problem when Q-balls form after the 
Affleck-Dine baryogenesis. Then, we study if the LH U flat direction can explain the baryon 
asymmetry in Sec. 5. The final section is devoted to summary and discussions. 

2 Moduli problem and Thermal inflation 

In this section, we review the moduli problem and the thermal inflation in gauge-mediated 
SUSY breaking models. 
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2.1 Moduli problem 


In gauge-mediation models, the moduli mass is much lighter than the electroweak scale, 
m v <C (9(100) GeV, since it acquires mass of the order of the gravitino mass through 
Planck-suppressed interactions. If the moduli mass is lighter than (9(100) MeV, its lifetime 
is larger than the present age of the universe. Therefore, the moduli energy density must 
be smaller than the dark matter density, /?dm — 0.26p cr , where pdm and p cr denote the 
present dark matter density and the critical density, respectively. The ratio of the critical 
density to the present entropy density is given by 


— ~ 3.6 x 10" 9 h 2 GeV, (1) 

so 

where So is the present entropy density and h denotes the present Hubble parameter 
in units of 100 km/sec Mpc 1 . Furthermore, the energy density of moduli with mass 
m v > 200 keV is more severely constrained from the observation of X-ray background 
spectra if the moduli held decays into photons via a coupling with the kinetic terms of 
the gauge fields [8, 9]. Hereafter, we focus on the case of m v < 200 keV so that the coherent 
oscillation of the moduli held can be the dark matter and explain the dark matter density. 

As mentioned in the introduction, the oscillating energy density of the moduli held 
easily exceeds the critical density. The VEV of the moduli held is expected to be of the 
order of the Planck scale during and after the primordial inflation since it is a singlet held 
without any symmetry enhanced points. When the Hubble parameter decreases to the 
moduli mass scale after the primordial inflation, the moduli held starts to oscillate and 
its energy density scales as a -3 , where a is the scale factor. In the case that the reheating 
after the primordial inflation occurs after the onset of the moduli oscillation, the ratio of 
its energy density to the entropy density is estimated as 


Py_ ^ 9rh / Vo \ 
Si ~ 8 V M p1/ 


~ 1.3 x 10 5 GeV 


T, 


10 6 Gev) 



( 2 ) 


where s t is the entropy density, Trh is the reheating temperature and M p \ (~ 2.4 x 
10 18 GeV) is the reduced Planck scale, r / 0 denotes the initial amplitude of the moduli 
oscillation, which is expected to be of the order of M p \. On the other hand, in the case 
that the reheating occurs before the onset of the moduli oscillation, the ratio is given by 


Py 


Si 


!/2 rnfal 

2tt 2 /45^T 3 


~ 1.2 x 10 6 GeV 


m v 

200 keV 



( 3 ) 


where g * is the effective number of degrees of freedom. T v is a temperature at the onset 
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supermultiplets Q u d L e H u H d X 


Z 4 charge 220020 2 1 

Table 1: Assignment of Z 4 charge. 


of the moduli oscillation and is defined as 



( 4 ) 


where we use g* = 229. The moduli-to-entropy density ratios are conserved until present 
unless entropy production occurs. From Eqs. (1), (2) and (3), one can see that the 
produced moduli energy density is much larger than the critical density. This is called 
“the cosmological moduli problem”. 

2.2 Model of thermal inflation with symmetry 

In this subsection, we explain a model of the thermal inflation with Z 4 symmetry and 
show that it is excluded because of a late-time decay of “i?-axions”. 

The thermal inflation [13] is an attractive solution to the cosmological moduli problem. 
It requires a singlet held X under the standard gauge symmetries, which is called “a 
haton”. We assume that it has a relatively hat potential by imposing a Z 4 symmetry 1 
and has interactions with particles in the thermal bath: 



( 5 ) 


where Ax and is dimensionless coupling constants and X is the supermultiplet of the 
haton held. £ and £ are massless gauge charged helds in the thermal bath. C is a 
constant term which cancels out the vacuum energy and is related to the gravitino mass 
by |C| ~ m 3 / 2 Mpj. We can ignore higher dimensional terms since the VEV of the haton 
is much smaller than the Planck scale. 

Motivated by solving the g problem in the Higgs sector, we assign Z 4 charge to the 
minimal SUSY Standard Model (MSSM) particles, as shown in Table. 1. We assume that 
the A-parity is also conserved in the MSSM sector. The Z 4 symmetry allows the following 
term [14, 15]: 



( 6 ) 


1 Even if we impose Z n (n > 4) symmetry to the superpotential, the following discussion is almost 
the same. 
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where H u and H d denote up-type and down-type Higgs supermultiplets. is a dimen¬ 
sionless coupling. Although the //-term is absent due to the Z 4 symmetry, the coupling 
of Eq. (6) generates an effective /i-term, which is necessary for the electroweak symmetry 
breaking, since the flaton acquires the VEV at low energy as we will explain below. 

Let us focus on dynamics of the flaton held in the early universe. At finite temperature, 
the baton potential is given by 2 

V(X) ~V 0 + ( cT 2 - mjS)|X| 2 + (^X* + h.c.'j + hh!|X|«, (7) 

where we use the same letter X to denote the complex scalar component of the baton 
supermutiplet. c is of the order of the square of the coupling constant g £. Here, we add the 
negative mass term (—ml < 0) induced from gauge-mediated SUSY breaking ebects [29], 3 
which is essential for the thermal inhation to work. When the baton is near the origin 
and cT 2 > m§, the baton held is trapped at the origin. 4 Then, the vacuum energy Vo 
becomes to dominate the energy density of the universe and exponential expansion of 
the universe, called “thermal inhation”, occurs. The thermal inhation lasts until the 
temperature decreases to the critical value of T c ~ At the temperature of T c , 

the baton starts to oscillate around the true minimum given by 

and the thermal inhation ends. The baton continues to oscillate until it decays, which 
releases huge entropy, and then the radiation dominated universe is realized. This is the 
mechanism of the thermal inhation, which dilutes the moduli relic by the huge entropy 
production of the baton decay. 

At low energy, the baton held X acquires the VEV and is decomposed as follows: 

- Y = ( M+ 7u) exp (y!i i)' (9) 

where \ and a are canonically normalized real scalar belds. The component \ corresponds 
to the baton, which starts to oscillate after the thermal inhation. If C was absent in 

2 Note that since particles which couple with the flaton acquire mass about g% ( X ), the flaton acquires 
the thermal mass term only around the origin of the potential, i.e., only when g j ( X) -C T. 

3 Although the flaton is a singlet field, the soft SUSY breaking mass term arises because it has a 
coupling with massless gauge charged fields, £ and £. 

4 The initial condition is naturally explained by the effects of supergravity. A positive Hubble induced 
mass term, which dominates over the soft mass term in the early universe, sets the flaton field at the 
origin. After that, it is trapped at the origin by the thermal mass term. 
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Eq. (5), X would have a U( 1)r symmetry and the field a would have a flat potential. 
We then refer to the scalar field a as “an i?-axion”, which obtains small mass from the 
nonzero constant term C (~ The masses of the flaton x and the A-axion a are 

calculated as 


m 2 — 4m 2 


0) 




( 10 ) 


respectively. Note that m x m a since the flaton obtains mass of the order of the SUSY 
breaking scale (m x — (9(10 2 -10 3 ) GeV) and m 0 m 3 / 2 . 

Let us consider the decay process of the flaton x- The flaton can directly decay into 
the lightest Higgs boson particles via the interaction of Eq. (6) when m x > 2m^, 5 where 
nth denotes the Higgs boson mass. The rate of this process is given by [16, 17] 


T(x ->■ 2h ) - 


1 

167T 


V M pi m x 



167r \ m x/ T7 2 ’ 


( 11 ) 


where /x = \X^\M 2 /M p \. Note that the flaton can decay into two higgsinos if it is kinemat¬ 
ically allowed, which may lead to overproduction of the LSP. We therefore assume that 
such a decay is kinematically forbidden, i.e., m x < 2/x. 

The flaton can also decay into i?-axions since m x > m a - At the tree level, the rate of 
this process is estimated as [16, 17] 


T(x ->• 2a) ~ 


1 m x 
647T M 2 


( 12 ) 


The produced i?-axions decay into Standard Model particles through one-loop diagrams, 
which is induced by the Yukawa coupling in Eq. (5) and the coupling of Eq. (6). In 
particular, i?-axions decay into photons and its rate is estimated as [16, 17] 


T(a —$■ 2y) ~ 


49 

72tt 



M 2 ’ 


(13) 


where a em is the fine-structure constant defined as a em = e 2 /47r. If loop running particles 
are charged under SU(3)c and m a > 1 GeV, the A-axion can also decay into gluons. 

As seen above, the flaton decays into both the Standard Model particles and R- 
axions when its decay rate becomes comparable to the Hubble parameter. The produced 
Standard Model particles thermalize immediately and reheating occurs. However, the 
produced i?-axions cannot thermalize soon since they have only one-loop suppressed in¬ 
teractions with particles in thermal bath. Later, R -axions decay into the Standard Model 

5 When m x < 2mh, the flaton decays into the Standard Model particles through one-loop diagrams. 
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particles and reheating occurs again. From Eq. (13), the reheating temperature when 
i?-axions decay, T F “ H , is estimated as 


T a ~ 
1 RH — 


90 

7T 2 g, 


1/4 


\J r(a — > 


“ 4.5keV|A,| 1/2 (^)( I ^) : 


( m 3 / 2 \ 3 / 4 

200 keV / 


(14) 


where we use g * = 3.36. It is found that the A-axions decay at the temperature much 
lower than (9(10) MeV. This implies that it destroys light elements synthesized at the 
BBN era and spoils the success of the BBN. 


2.3 Model of thermal inflation with a linear term 


One way to avoid the late-time decay of A-axions in the Z 4 symmetric model is a prohi¬ 
bition of the flaton decay into A-axions, which can be realized by adding a Z 4 symmetry 
breaking term [16, 17], 

SW = aX, (15) 


into the superpotential. Note that a is a parameter with mass dimension 2. As we will 
explain below, this term leads to degeneracy of masses of the flaton and the A-axion, and 
the flaton decay into i?-axions is forbidden kinematically. 

Here, we should comment on another motivation to introduce the Z 4 symmetry break¬ 
ing term. When the superpotential has the Z 4 symmetry as we consider in the previous 
subsection, the flaton has four degenerate minima and the model is faced with the do¬ 
main wall problem after the thermal inflation. However, the Z.\ symmetry breaking term 
of Eq. (15) resolves the degeneracy of the vacua, which leads to fast collapse of domain 
walls. The problem can be solved if they collapse before their energy density dominates 
that of the universe. Quantitatively, the condition to solve the domain wall problem is 
given by [9, 15, 30] 


m 3/2 m x M 

Ml 


(16) 


Hereafter, we consider the case that |cx| is large enough to satisfy Eq. (16). 
The flaton potential is now given by 


V{X) = Ho 




mm 


a* A 


X V 3 


M, 


X 6 + h.c. 


pi 


1 CX 


X V 4 


4M 3 


X 4 + h.c. + 


|A 


x 


K 


■\x[ 


(17) 












The negative mass term is induced from the gauge-mediated SUSY breaking effect and is 
applicable only for the the flaton VEV smaller than the messenger mass scale. When the 
flaton VEV becomes larger than the messenger scale, the negative soft SUSY breaking 
term becomes suppressed at large amplitudes. In the following, we focus on the parameter 
region where the flaton VEV at the true minimum is larger than the messenger mass scale 
and the negative mass term is negligible. In order to obtain the effective //-term of the 
order of the electroweak scale via the coupling of Eq. ( 6 ), the flaton VEV should be larger 
than 0(1O 10 ) GeV. In this parameter region, the trilinear term dominates over both the 
linear term and the quartic term at the true minimum. The flaton VEV at present and 
the vacuum energy at the origin are thus estimated as 

{x) = M~ , V - H 2 , (18) 

where Vq is determined by requiring that the vacuum energy vanishes at the true minimum. 
The masses of the flaton and the Y-axion are now almost the same and are given by 

o f\ a \ 2 \^x\\ 1/3 

ro * m * HhWpTJ ' (19) 

Note that the flaton decay into Y-axions is forbidden since their masses are degenerate. 

Taking the thermal effects into account, the flaton potential around the origin is 
rewritten by 


V{X) ~ (cT 2 -ml) 


X — 


2 a*C 


(cT 2 - m 2 )Ml 


+ U 0 


4\a\ 2 C 2 


(cT 2 -m 2 )M^ 


and the local minimum before and during the thermal inflation is given by 

2\a\m 3/2 




cT 2 - ml' 


( 20 ) 


( 21 ) 


Then, one can find that the flaton does not sit at the origin due to the Z 4 symmetry 
breaking term a. This implies that the thermal inflation ends when the VEV of X 
becomes so large that it decouples from the thermal plasma, i.e., when the temperature 
decreases to T c ~ c 1//2 (|Y|). Although the vacuum energy during the thermal inflation 
also deviates from Vo, the deviation is estimated as 5V/V ~ my 2 /cT 2 and is negligible. 

We turn to estimate how much entropy is released by the flaton decay after the thermal 
inflation. We define a dilution factor as the ratio of the initial to the final comoving entropy 
density. It is calculated as [16, 17] 


s f a f 4 

A = A-4 = 1 + 


U n 


S i a i 


32*745 


( 22 ) 


RH 
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where Sf represents the entropy density after the thermal inflation and T£ h denotes the 
reheating temperature when the flaton decays. T^ H is calculated from its decay rate and 
has to be higher than 0(10) MeV. 

From the dilution factor, one can estimate the moduli density parameter. Hereafter, 
we refer to the moduli produced before the thermal inflation as “Big-Bang moduli”. Using 
Eq. (2), the density parameter of the Big-Bang moduli is given by 




1 Pri, BB So h 2 


A s. 
~ 0.026 


Pci 

T, 


RH 


Un 


10 6 GeV / \ 10 24 GeV 4 


rpX 

1 RH 

10 MeV 


T r 


1 TeV 


Vo 


M 


pi 


■ (23) 


where the subscripts “BB” represent the Big-Bang moduli. If the reheating of the primor¬ 
dial inflation occurs before the moduli oscillation, one can obtain the density parameter 
using Eq. (3). It is found that the abundance of the Big-Bang moduli density can be 
sufficiently diluted by the thermal inflation. 

Unfortunately, the thermal inflation gives a secondary source of the moduli. We call 
the secondary produced moduli as “thermal inflation moduli”. The potential of the moduli 
held during the thermal inflation is given by 


V(IJ) = 




1 


+ y h ' 2 (v~Vo) 2 


(m + c h H ) 77 - 


c h H 2 


717.2 + C H H 2 


VO + 


(24) 


where ch is a coefficient of the Hubble induced mass term, which comes from the finite 
vacuum energy during the thermal inflation. 170 is a minimum determined by the Hubble 
induced mass term and is expected to be of the order of the Planck scale. Equation (24) 
clarifies that the moduli held sits at the minimum depending on the Hubble induced mass 
term during the thermal inflation. After that, it starts to roll down to the true minimum 
with the amplitude of 77 ~ Ch {Hxh/m^) 2 770 , where H th is the Hubble parameter at the 
end of the thermal inflation. The ratio of the moduli to the entropy density at present is 
therefore given by 

P’i.TH _ faTH _ <).V| ,25-, 

Sf 4V37?„ 24’ 

The density parameter of the thermal inflation moduli is estimated as 




Pr ),TH SqH 2 

Sf Pci 


~ 0.484 


Vo \ 
10 24 GeVV 


( \- 2 / 4H A 

V 200 keV / V 10MeV / 



(26) 
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Note that the thermal inflation moduli is produced more abundantly as the vacuum energy 
Vo becomes higher. This is contrary to the Big-Bang moduli. 

From Eqs. (23) and (26), it is found that D^bb^ 2 + ^,th^ 2 0(0.1) for Vo ~ 

0(1O 24 ) GeV 4 , T c ~ 0(1) TeV and T^ H ~ 0(10) MeV, and that the thermal inflation can 
solve the moduli problem. The relic of the moduli could be the cold dark matter. 

3 Affleck-Dine baryogengesis and Q-ball 

As mentioned in the introduction, the thermal inflation dilutes not only the moduli den¬ 
sity but also the baryon asymmetry which may be produced beforehand. Therefore, some 
viable mechanisms to produce sufficiently large amount of baryon asymmetry are needed 
for a cosmological consistent scenario. The Affleck-Dine baryogenesis [20, 21] is a promis¬ 
ing candidate in the framework of SUSY. In the Affleck-Dine mechanism, flat directions 
rotating in the complex plane produce baryon number and they decay into the Standard 
Model particles in the early universe. In gauge-mediated SUSY breaking models, however, 
the Affleck-Dine field forms non-topological solitons called “Q-balls” and the produced 
baryon charge is confined into Q-balls [22, 23, 24, 25]. We briefly review the Affleck-Dine 
baryogenesis and Q-ball properties in this section. 

3.1 Affleck-Dine baryogenesis 

Supersymmetry predicts a lot of flat directions which have vanishing potentials in renor- 
malizable levels. In particular, the MSSM contains flat directions carrying baryon (and/or 
lepton) charge. They are referred to as “Affleck-Dine (AD) fields”. One example of them 
is the udd flat direction. The udd flat direction is a linear combination of right-handed 
squark fields with the same amplitude, for example, uf = df = and df = ^<3>, 
where u and d are the up-type and down-type right-handed squarks, respectively. The 
superscripts and the subscripts show color and family indices, respectively. In this direc¬ 
tion, D-terrn potentials indeed vanish and it also has no renormalizable terms in F-terms 
if the A-parity is conserved. In the following, we explain the Affleck-Dine mechanism 
without restricting to the udd flat direction. 

The flat directions acquire the potential from the SUSY breaking effects and the scalar 
potential is expressed as 

U SB ($) = U gauge + U grav = Mp (log + m 2 3/2 A + K log |<f>| 2 , (27) 

\ mess / \ ^ / 


li 




where <h denotes the AD held. The first term comes from gauge mediation effects [6, 23] 
and is applicable only for |<3>| M mess , 6 where M mess is the messenger mass scale. Mp is 
related to SUSY breaking scale as follows: 


Mp — 


9 —Vw 

47r 


(28) 


where g generically represents gauge couplings of the Standard Model and F denotes the 
SUSY breaking F-term. The parameter k is determined from couplings between the SUSY 
breaking sector and the messenger sector and satisfies k < 1. The SUSY breaking scale 
is constrained from below by the relatively heavy Higgs boson mass at around 126 GeV, 
which acquires the radiative corrections from the stop mass. Since the masses of scalar 
particles are proportional to the parameter A = kF/M messi we obtain a lower bound of A 
as follows [31, 32]: 


A = 


kF 


M n 


> 5 x 10 5 GeV. 


(29) 


In combination with M^ ess > kF which is necessary to make masses of messenger scalars 
positive, \JkF > A > 5 x 10 5 GeV must be realized. Therefore, Mp is restricted from 
below, 

M F > 4 x 10V /2 GeV. (30) 


Since the SUSY breaking scale F is related with the gravitino mass as F ~ v / 3^3/2Af p i, 
Mp is restricted from the above as follows: 


Mp< 2 x 10 6 GeVg 1/2 ( m3/2 ) ^ . (31) 

y V 200 keV / V ’ 

The second term of Eq. (27) comes from the gravitational mediation effects including 
one-loop corrections [33]. M* is a renormalization scale. The parameter K comes from 
1-loop effects and its absolute value is typically in the range of 0.01 to 0.1. If gaugino 
contribution to 1-loop effects is larger than that of the Yukawa interactions, the sign of 
K is negative, and vice versa. Notice that the gravitational mediation effects dominate 
over the gauge mediation effects above a threshold value (p eq which is given by 


0 eq ~ 3.2 x 10 14 GeV 


/ "U/2 X" 1 
\200 keV / 


( M f V 

\4 x 10 4 GeV ) ’ 


(32) 


for M mess ~ 5 x 10 5 GeV. 

6 For |<F| <C M m ess, the first term in Eq. (27) is replaced by a soft SUSY breaking mass term, 

to 1 usyI^| 2 i where msuSY ~ D(0.1-l) TeV. On the other hand, it becomes logarithmic for |4>| M mess 

since a large amplitude of the AD field leads to suppression of the transmission of SUSY breaking 

effects [6]. 
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There also exists thermal effects on the effective potential. When /1$| <C T ,where / 
is a coupling constant between thermal particles and the AD field, the thermal mass term 
arises and is given by 


A 1 ’ ~ / 2 r' 2 |<s | 2 


(33) 


A thermal effect via 2-loop diagrams also contributes to the effective potential for /|<f>| 
T as follows: 


K 


( 2 ) 


th 


o^T 4 


log 


|$| 2 

J>2 


(34) 


where a g = g 2 /47r [34, 35]. 

The flat directions are also lifted by non-renormalizable superpotential and we gener- 
ically write it as 


Wnr — 


nM p "f 3 ’ 


(35) 


where A$ is a coupling constant and n (> 4) is an integer. For example, in the case of 
udd flat direction, n = 6. The superpotential of Eq. (35) leads to the non-renormalizable 
terms in the AD field potential as 


l / NR( < f ) ) — \<s>a\m3/2 


cj> n 


I (jjpn—2 

+ h-C. + VP ^2n-6 • 


(36) 


Here, we introduce a dimensionless parameter a\ of Oil), which depends on higher dimen¬ 
sional Kaliler potentials. The terms proportional to a\ are induced from the gravitational 
effects and called “A-terms”. 

In addition, the AD field acquires Hubble induced terms in the early universe. The 
vacuum energy causing the expansion of the universe largely breaks the supersymmetry 
and lifts the potential of the flat direction. Since the vacuum energy is related to the 
Hubble parameter, these additional terms are given by 


V H 


— chH~\&\ 2 T A<j mhH 




+ h.c., 


(37) 


where Ch and an are dimensionless parameters of 0(1). Hereafter, we consider the case 
of ch > 0 to realize the AfHeck-Dine baryogenesis. 

Next, let us explain dynamics of the AD field in the early universe. Its potential 
is given by Eqs. (27), (36) and (37). During the primordial inflation, the AD field sits 
down at the minimum determined by the negative Hubble induced mass term and the 
non-renormalizable terms in Eq. (36). Its VEV is estimated as 


4>| ~ (HM«- VIAtl) 17 '” 2) 


(38) 
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and the phase of the AD held is determined by the Hubble induce A-terrn. After the 
inflation end, the energy density of the universe is dominated by the oscillating inhaton 
and the Hubble parameter decreases. When H ~ W($)/$, the AD held starts to roll 
down to the origin of the potential. At this time, the phase of the AD held also starts 
to rotate in the complex plane producing baryon number, Ub — Im[$*$]. The produced 
baryon number density at the onset of the oscillation of the AD held is roughly estimated 
as 


n B I ad 


if -1 Im 


dV 


$ 


AD osc 


^3/2 I'M |^osc 
hToscMpV 3 


( 9 ^) 14 

\ -^osc / 


osc ^ osc 


= eH 1$ r 

— c -*-* OSC I ^OSC I 1 


(39) 


where we use Eq. (38) and dehne the ellipticity parameter e = 1713 / 2 /H osc . As the oscilla¬ 
tion amplitude of the AD held decreases due to the Hubble friction, higher order terms 
become irrelevant. Therefore, the baryon number violating operators, the A-terms, pro¬ 
duce baryon number only at the onset of the oscillation. After that, the comoving density 
of the baryon number is conserved. The produced baryon number is finally converted 
into the Standard Model particles by the decay of the AD held and the present baryon 
asymmetry can be explained. This is the mechanism of the Afheck-Dine baryogenesis. 


3.2 Q-ball 

It is known that the AD held oscillating in the potential of Eq. (27) feels spatial instabilities 
and forms into non-topological solitons, Q-balls [22]. Let us explain a condition that Q- 
balls are formed after the Afheck-Dine baryognesis. The configuration of the Q-ball is 
determined by the condition of minimizing the energy with conserved baryon charge, 
where the energy and the baryon charge are given by 


E 


d 3 x(\<5>\ 2 + |V<L| 2 + H(|<L|)) , Q = 2 


d 3 x Im 




(40) 


respectively. Here, we assume that the AD held carries an unit of baryon charge for 
simplicity. The scalar held configuration is obtained by minimizing 


E u = E + Ld 



d 3 x Im 



(41) 


where oj is a Lagrange multiplier. E u is rewritten by 


E U1 = / d 3 x 


$-fca<h| i -a; 2 |$| 2 +|V$| 2 + H(|$|) +wQ. 


(42) 
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The time dependence of $ is determined as <E>(x, t) = p(r)e lut / \/2 from the first term 
in order to minimize E w . We assume that the stable field configuration is spherically 
symmetric, which leads to the following equation of the field configuration: 


d 2 


2 d 


5^ + ^ + ^ 


|v^) = c. 


(43) 


The boundary condition is ^(0) = 0 and <^( 00 ) = 0 in order to obtain a smooth and local 
configuration. The solution with the boundary condition exits for 


u 2 0 


mm 


' W(<P) 

L 


< 


1 / 3 = 1 / 30^0 


d 2 V{0) 

dip 2 


(44) 


This inequality requires that there exists a field value where the potential is flatter than 
the quadratic potential. This condition is satisfied if the potential of the AD field is 
dominated by the first term of Eq. (27). Even if the second term dominates over the first 
one, the condition of Eq. (44) is satisfied for the case of K < 0. Numerical simulations 
have shown that AD fields which oscillate in the potential satisfying Eq. (44) actually 
form into Q-balls [25, 26, 27]. ft has also been shown that (almost) all baryon number is 
confined into Q-balls. 

The profile of the Q-ball depends on the potential of the AD field. When the gauge- 
mediation effect dominates over the gravity-mediation one, formed Q-balls are referred to 
as “gauge-mediation type Q-balls”. The field configuration of the gauge-mediation type 
Q-ball is determined by solving Eq. (43) and is approximately given by [36] 


giujt 

H r ) - -j= x 

where uj and p 0 are given by 


for r <R = vr/cu 
0 for r > R 


(45) 


u — V2nM F Q 1 / 4 , ip 0 ~ M f Q 1 / 4 . 


(46) 


The energy of the Q-ball is calculated from Eq. (40) and is estimated as 

E ^ (47) 

o 

One can find that the Q-ball energy per unit charge (~ dE/dQ ) is smaller for larger 
Q. When it is smaller than the proton mass, i.e., dE/dQ rs./ MpQ- 1 / 4 < 1 GeV, Q-balls 
cannot decay into nucleons. The charge of the Q-ball can be determined by numerical 
simulations and is given by [27] 

(48) 
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where /? is a numerical coefficient and is determined as /3 ~ 6 x 10” 4 .' Hereafter, we use 
4> as the amplitude of the AD held and 0 OSC denotes the amplitude at the onset of the 
oscillation, |<f> osc |. Note that 0 OSC should be smaller than 0 eq to use Eq. (48). 

In the case of </> osc > </> eq , the gravity mediation effect dominates over the gauge 
mediation one and one can consider two types of scenarios: K > 0 and K < 0. When 
K < 0, the second term of Eq. (27) satisfies the condition for Q-ball formation. Q-balls 
formed in this potential are referred to as “new type Q-balls” [26]. The held configuration 
is approximately given by a Gaussian function [33]: 




-r 2 /2R 2 iut 


(49) 


where R, uj and <p Q are given by 


R ~ 


AT / 2 


m 3/ 2 


, u ~ m 3 / 2 , <^o — 


® ) 3H m 3/2 Q^ 


(50) 


respectively. The energy of the Q-ball is given by 


E ~ m 3/2 Q. 


(51) 


This type of the Q-balls is stable against the decay into nucleons since dE/dQ ~ m 3 /2 < 
1 GeV in gauge-mediated SUSY breaking models. The charge of the new type Q-ball is 
given by [37, 38] 

Q~p(—') (52) 

\m 3 / 2 J 

where (3 ~ 2 x 10 -2 . When K > 0, on the other hand, the condition of Eq. (44) is 
not satished and Q-balls are not formed. In this case, the oscillation of the AD held 
remains homogeneous and its amplitude decreases as </> oc a -3 / 2 after it starts to oscillate. 
However, Q-balls are formed when its amplitude decreases to <T q and the potential of 
the AD held becomes dominated by the gauge mediation effect. This type of Q-balls is 
referred to as “delayed type Q-balls” [27]. The prohle and properties of the Q-ball are 
the same as those of the gauge-mediation type Q-ball (see Eqs. (45), (46) and (47)), while 
the charge of the delayed type Q-ball is given by 

(53) 

7 Precisely speaking, the numerical coefficient /3 depends on the orbit of the AD field in the complex 
plane. /3 ~ 6 x 10~ 4 for a circular orbit (e = 1), while /? ~ 6 x 10 -5 for an oblate orbit (e < 0.1). 
Hereafter, we use f3 ~ 6 x 10~ 4 for simplicity since our results do not change significantly. 
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The thermal logarithmic potential of Eq. (34) also satisfies the condition for Q-ball 
formation. If the AD field starts to oscillate by the thermal logarithmic potential, one 
can obtain the charge of the Q-ball by replacing Mp with T* in Eq. (48), where T* is the 
temperature at Q-ball formation. The profile and properties become the same as those of 
the gauge-mediation type Q-balls when the temperature decreases sufficiently. 


4 Moduli problem and Baryon asymmetry with Q- 
ball formation 

In gauge-mediated SUSY breaking models, Q-balls are formed during the oscillation of 
the AD field. Almost all baryon number is absorbed into Q-balls even if huge baryon 
asymmetry is produced by the Affleck-Dine baryogenesis. If Q-balls are stable, baryon 
charge is released from Q-balls via evaporation [39]. In the case of unstable Q-balls, 
the decay of Q-balls can release baryon charge. However, it must occur before the BBN 
epoch in order to explain the baryon asymmetry. In this section, we show that the baryon 
number production is incompatible with the moduli dilution in both cases updating the 
analysis in [28]. 


4.1 Stable Q-ball 


First, we consider stable Q-balls and briefly review that it is difficult to explain the 
baryon asymmetry in this case. The Q-ball is stable against the decay into nucleons when 
dE/dQ < 1 GeV. At finite temperature, however, baryon charge can evaporate from Q- 
ball surface since the free energy of the AD field in the thermal plasma is smaller than 
that of Q-balls. The total evaporated charge of gauge-mediation type Q-balls is given 


by [27, 40] 


A Q ~ 10 16 


f m susy 
V ITeV 


- 2/3 


M P 


4 x 10 4 GeV 


- 1/3 


Q l,1 \ 


(54) 


where mgusy denotes a soft SUSY breaking mass scale. In the case of new type Q-balls, 
it is given by [26, 40] 


A Q ~ 10 20 


\K\_\ V3 / ^ 3/2 
0.01 ) \ 200 keV 


-V 3 / msusv 
V ITeV 



(55) 


Since A Q/Q becomes smaller as Q increases, it is more difficult to extract baryon charge 
from Q-balls with larger baryon number. In order to produce huge baryon number, </> osc 
should be large, which increases the baryon charge of the Q-ball. Then, the released 
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baryon number generally becomes smaller as (j) osc increases. In the case of the delayed 
type Q-ball, however, the baryon charge of Q-balls does not increase even though </> osc 
becomes much larger than 0 eq . Hence, the delayed type Q-ball seems to be able to most 
effectively provide baryon charge outside Q-balls. We thus focus on the delayed type 
Q-ball in the following. For the delayed type Q-ball, A Q/Q is estimated as 


^^4 iq —18 f^SUSY\~ 2/3 ( M f \ 4 / m 3/ 2 yi/3 

Q ~ V 1 TeV / V 4xl ° 4GeV / v200 keV/ 


(56) 


where we use Eq. (32). 

The finally provided baryon asymmetry is calculated from 


n B _ n B A Q 
s s Q 


n B 

Pr], bb 


77 osc 


Pr), BB 
S 


A Q 
0 Q ’ 


(57) 


where n B /s is the ratio estimated without considering Q-ball formation. When </> osc > 0 eq , 
the ratio of the baryon number to the Big-Bang moduli energy density is estimated as 


n B 

n B 

3 M 2 # 2 


( </wV 

( m y 2 

Pr), BB 

3 M 2 TT 2 

77 OSC pi 

AD osc P "- BB 

H 2 

77 osc osc 

[m p J 

{ MpJ 


2 / 0o,o y / vo_\ 2 

m 3 /2 V M pl/ \ V/p| / 


(58) 

(59) 


where we assume H osc ~ rri 3/2 in the second line, which is the case of new or delayed 
type Q-balls. Here, we assume that the total energy density of the universe scales as a~ 3 
between the eras when the AD field starts to oscillate and when the moduli held starts to 
oscillate, since the moduli held starts to oscillate before the reheating. 8 From Eqs. (56) 
and (59), the baryon asymmetry is expressed as 


Y b 


rs./ 


2 x 10- 23 ( 

v 0.12 

f rasusv W 2 / 3 

V 1 TeV ) 


M p \ 

M f 


-2 


4 x 10 4 GeV 


/ m 3 / 2 \ 8 /3 
\ 200 keV / 



(60) 


where we use the Big-Bang moduli density after the entropy dilution due to the thermal 
inflation and assume that the abundance of the thermal inflation moduli is negligible. 

8 Fgauge > F t ( h 2) and F gauge > F grav lead to 4> osc /M pl > (T RH /M F ) 2 and (j> OS c/M p \ < Mj./m 3/2 M p i. 
Combining these relations, Trh < M F / sjrri 3 / 2 M p \ should be realized. Since Mp < m 3 / 2 M p \ , the re¬ 
heating temperature is estimated as Trh Tj, which implies the reheating occurs after the moduli 
oscillation. 
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Here, we assume that the thermal effects on the effective potential is negligible. The 
stability condition, dE/dQ < 1 GeV, leads to the upper bound on the gravitino mass, 
777 . 3/2 < 1.4 GeV. One can find that the estimated baryon asymmetry is too small to 
explain the present one, Y F — 8 x 10 11 , even if 0 OSC ~ M p \ and M F taken as the 
lower bound of Eq. (30), which puts the more severe upper bound of the produced baryon 
number than [28]. Other scenarios, including ones where Q-balls are formed by the thermal 
logarithmic potential, cannot explain the present baryon asymmetry, either [28]. 


4.2 Unstable Q-ball 


We show that the situation does not become improved even in the case of the unstable 
Q-ball. It can decay into nucleons when dE/dQ > 1 GeV and release baryon charge from 
its surface. The emission rate is determined by the Pauli blocking effect on its surface 
and is given by [41, 42, 43] 


dQ ~ (2 u) 3 A 
dt 967T 2 


(61) 


where A is the surface area of the Q-ball. In order to avoid destroying light elements due 
to the decay products of Q-balls, the decay temperature is constrained as follows: 


Tde c 



Mpi > 0(10) MeV. 


(62) 


In gauge-mediated SUSY breaking models, unstable Q-balls correspond to “gauge- 
mediation type Q-balls”, where the AD held begins to oscillate by Ugauge- In this case, of 
course, Uauge > Vgrav should be satisfied, which leads to 0 OSC < 0 eq - Moreover, dE/dQ > 
1 GeV leads to 


— < 1.2 x 10 
M pl ~ 


-5 


/ m f y 

VlO 6 Gevy 


(63) 


From Eq. (62), the amplitude of the AD held at its oscillation is constrained as 


M 


< 1.7 x 10" 6 


pi 


/ M f \ 6/5 

VlO 6 Gevy 


The resulting baryon asymmetry is then estimated as 


(64) 


Y b = 




kb Pr ),bb 


Pri, BB S 

1 0 f / VO \ 2 / ^ 3/2 \ / M F 

' V 0.12 y \W p ij \ 200 keV/ V ^ 6 GeV 


-4 


0 osc 

Mpi 


4 


1 


(65) 
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where we use Eq. (58) and H osc ~ Mp/(p osc . Since (j) osc is restricted from above by Eq. (64), 
the baryon asymmetry has an upper bound of 


Y r < 8.5 x lO” 24 


BB ^ 2 

0.12 


Vo 


M, 


pi 


-2 


m 3 /2 \ 


Me 


1 


4/5 


200 keV/ V 10 6 GeV / 


( 66 ) 


Note that Eq. (63) and (f> osc < 0 eq are satished when M F > 9.0 x 10 4 GeV and m 3 / 2 < 
0(200) keV. One can find that the estimated baryon number is too small to explain the 
present baryon asymmetry. This results from the smallness of the amplitude of the AD 
held to make Q-balls unstable. 

In the above discussion, we neglect the thermal logarithmic potential. Next, we assume 
that the reheating completes before the moduli oscillation and that the AD held begins 
to oscillate by the thermal logarithmic potential before the reheating. 9 Even in this case, 
Eq. (58) is applicable when 0 OSC is near the Planck scale since the AD held immediately 
dominates the total energy density after the reheating. The charge of the Q-ball is given 
by Q ~ /d(0osc/Tosc) 4 , where T osc is the temperature at onset of the AD held oscillation 
and is estimated as T osc ~ (MpiT RH i/ osc ) 1//4 5 with H osc ~ T 2 SC /0 OSC . dE/dQ > 1 GeV leads 
to 


M 


< 3.7 


pi 




10 6 GeV 


2/3 


T-i 


RH 


\ 2/3 

6 x 10 11 GeV ) 


From Eq. (62), the amplitude of the AD held at its oscillation is constrained as 


M t 


< 1.0 


pi 


M F \ 
10 6 GeV ) 


2/15 


\ 2/3 

6 x 10 11 GeV ) 


iRH 


(67) 


( 68 ) 


The resulting baryon asymmetry is estimated as 


Y b = 


n B Prj, BB 
Pr), BB S 

~ 2.9 x 10” 50 


h^»7,BBh 2 

0.12 


Vo 


M, 


This has an upper bound of 10 


Yb < 8.5 x 10” 24 




,h 2 


0.12 


pi 


Vo 


-2 


RH 


m 3/ 2 \ 

200 keV / VM, 


T, 


i p i 


-4 


M, 


M, 


pi 


-2 


m 3 /2 \ 




1 


pi 


4/5 


200 keV / V 10 6 GeV / 


(69) 


(70) 


where we use Eq. (68). Note that Eq. (67) is satished when M F > 9.0 x 10 4 GeV. The 
resulting baryon asymmetry is too small to explain the observed baryon asymmetry. This 
results from the early oscillation of the AD held by the thermal logarithmic potential. 


9 Vh ~ Egauge and Vh ^ Egrav lead to (j) osc /M v \ < (T rh /Mf) 2 and A>sc/T/pi < (Tr H /to 3 / 2 M p i) 1 / 2 . 
These are satisfied if the reheating completes before the moduli oscillation. 

10 By taking into account of the ellipticity parameter e ~ m 3 / 2 /H osc , we obtain more severe upper 

bounds of Eq. (66) and Eq. (70) than the estimated value in [28]. 
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5 LH U flat direction 


In the previous section, we focus on flat directions which have only soft SUSY breaking 
terms and non-renormalizable terms. However, there exists an exception. The LH U flat 
direction has the supersymmetric /i-term of Higgs supermultiplets in addition to the SUSY 
breaking terms. Moreover, one-loop correction gives a positive contribution (A" >0). This 
implies that the LH U flat direction does not form into the Q-ball. 

First, let us consider the case that the /i-term exists before the thermal inflation and 
show that it is still difficult to explain the baryon asymmetry. Assuming that the AD 
held begins to oscillate because of the /i-term before the thermal inflation, the ratio of 
lepton number to Big-Bang moduli density is given by 


n L 

2m 3/2 

( d'oscV / 

Pv, BB 

P 2 

\ 



-2 


J 


(71) 


where /i is the electroweak scale. Here, we use Eq. (58) which is applicable when the total 
energy of the universe scales as a~ 3 between the eras when the AD held starts to oscillate 
and when the moduli held starts to oscillate. Even if the reheating completes before the 
onset of the moduli oscillation, the AD held immediately dominates the total energy of 
the universe when 0 OSC ~ M p \ and the above estimation is valid in that case. The lepton 
asymmetry is partially converted into baryon asymmetry through the sphaleron process. 
The relation between the lepton and baryon number is given by 


n B __8 n L 

s 23 s 


The produced baryon asymmetry is estimated as 


(72) 


Y b ~ 6.1 x 1(T 20 


V 0.12 ) 


VO \ 2 ( ^3/2 \ / \l \~ 2 /0osc\ 2 

M p \) \ 200 keV / VlTeW \M P J ’ 


(73) 


and the resulting asymmetry is too small to explain the present asymmetry. This results 
from the early oscillation of the AD held by the /i-term whose scale is much larger than 
the gravitino mass scale, /i m 3 / 2 . 

Hereafter, we consider an alternative scenario where the /i-term is negligible at the 
onset of the oscillation of the AD held and is generated at the end of the thermal inflation. 
The coupling between the baton and Higgs supermultiplets (see Eq. (6)) prohibits the id- 
term before the thermal inflation and plays the role of generating the /i-term by the VEV 
of the baton after the thermal inhation. In this case, the scenario changes as follows. The 
LH U hat direction produces lepton asymmetry by soft SUSY breaking terms and forms 
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into Q-balls before the thermal inflation. After the thermal inflation, the generated /!- 
term violates the condition for the existence of the Q-ball. Q-balls decay and the absorbed 
lepton number is released. In the following, we show if the released lepton number can 
explain the observed baryon asymmetry. 

In order to generate the /i-terrn, we assume that the flaton X couples with the Higgs 
supermultiplets as Eq. (6). Then, the following F-terrn potential arises: 

= (74) 


where $ is the up-type Higgs scalar held, namely the AD held. For convenience, we 
introduce the p parameter defined as Eq. (74). Then, the baton potential is expressed as 


V(X) ~ Ho + 


a*X 


x 


Mr 


X 3 + h.c. + 


pi 


|A M | 2 |$| 

K 


W + 


|A 


x\ 


Ki 


■w 


(75) 


where we neglect the linear terms and the quartic terms in Eq. (17) (see the discussion 
in Sec. 2.3). Before the thermal inflation, the VEV of X is so small that /i-term is 
negligible for the dynamics of the AD held compared with SUSY breaking terms, which 
implies that the AD held forms into Q-balls. Then, the baton acquires the VEV after 
the thermal inhation, which gives the p-term to the AD held potential. If it breaks the 
condition for the existence of the Q-ball, the Q-ball decays. The lepton charge is then 
released to the thermal plasma. 

The /i parameter must increase from the outside to the inside of the Q-ball in order to 
violate the condition for the existence of the Q-ball. Namely, the following relation must 
be realized: 


a» 2 (v(|«„»,|)) ( ) 

<hin and $ out show the AD held values inside and outside the Q-ball, respectively. From 
Eq. (74), one can find that the /i parameter explicitly depends on |Y|. |V| is determined 
by the potential of Eq. (75), which depends on |$| via the interaction of Eq. (74). Thus, 
one can find that the p parameter depends on the AD held value. 

Since |$| ^ 0 inside the Q-ball, the coupling term between the baton and the AD held 
lifts the baton potential and \X\ inside the Q-ball is smaller than that outside the Q-ball. 
Hence, Eq. (76) is not satished at the tree level. Taking one-loop corrections into account, 
however, we find that the p-term can be steeper than a quadratic mass term for a larger 
VEV of the AD held |$|. It is expressed as 


V 


|A/M 4 

K 


^1 + K log 




(77) 
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where K > 0 and \K\ ~ 0(0.1-0.01). The parameter M is the VEV of the flaton. For 
different amplitudes of the AD held (|<E>+| > |<E>_|), the ratio of the //-term is estimated as 


p 2 (|$ + |) M% 
// 2 ( |$_|) Ml 


1 + K log ■ 




(78) 


where M + and M_ are the baton held values at |<E> + | and |$_|, and satisfy M + < M_. 
Precisely speaking, when the second term in the parenthesis is 0(1), the perturbation 
breaks down and one should solve renormalization group equations. For simplicity, we 
assume that the one-loop correction factor is 2 at most 11 and require 


1 ML 

2 


< i 


(79) 


to realize the condition of Eq. (76). Here, M in and M out correspond to the baton VEVs 
inside and outside the Q-ball. 

In order for the baton VEV not to be highly damped inside the Q-ball, the quartic 
term of the baton in Eq. (75) should be small enough. Thus, the AD held value in the 
Q-ball, (j) o, which is proportional to the oscillation amplitude in the case of the gauge- 
mediation type Q-ball, should be small. On the other hand, in order to produce huge 
baryon number, 0 OSC should be large. Hence, we focus on the delayed-type Q-ball because 
0o is independent of 0 OSC . 0o is determined by the held value, where the gauge mediation 
ebect is comparable to the gravity mediation ebect, and is estimated as 

2 


~ -^=/3 1//4 0 e q — 3.6 x 10 13 GeV ^ 


m 3 / 2 \ — 1 
200 keV / 


Mp 


4 x 10 4 GeV 


(80) 


Substituting this estimated value into the baton potential of Eq. (75), the baton VEV 
inside the Q-ball is obtained by solving the following equation: 


MIAyI 


M 


M 2 + 2 

in 1 " 


pi 


\K\^l 

K 


Mi + = 0. 


(81) 


In order to estimate the parameter a, we introduce a dimensionless parameter ( as follows: 


|ck| = (- 


M 


(82) 


pi 


Note that ( = 1 corresponds to the case without the coupling of Eq. (74) and that ( > 1 
is satisbed. Using Eqs. (18) and (82), one can bnd that the ratio of the baton VEV inside 
to outside the Q-ball is given by 

M, 1 

(83) 


M, 


out 


c 1/3 ’ 


Even if we assume that the correction factor is larger than 2, our result does not change significantly. 
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and the condition of Eq. (79) leads to 1 < ( < 1.7. By solving the equation of Eq. (81), 
one can obtain 

2|-V| 2 0o 


M = 

m 3(C-1)|A X | 2 

in terms of (. Then, from Eq. (82), the parameter a can be estimated as 


|a| = C 


3/2 


VI Vo 


(84) 


(85) 


L3(C — 1)J |A x | 2 M pl 

By using the expression of the fi- term as Eq. (74) and 1 < ( < 1.7, one can obtain the 
following constraint on the parameter cm 

1 1/2 

MI*out|)<fo 

MI*out|) 


a 


= (V3 


L3(C -1) 

> 1.3 x 10 16 GeV 2 


300 GeV / V 3.6 x 10 13 GeV 


( 86 ) 


Hence, when the symmetry breaking parameter a satisfies the above constraint, the con¬ 
dition of Eq. (76) can be realized including one-loop corrections. Hereafter, we assume 
the Z 4 symmetry breaking term of the order of |a| ~ 1.3 x 10 16 GeV 2 . 

From Eq. (18), the energy of the thermal inflation is estimated as 12 

/ mi \ 2 

Vn ~ 1.7 x 10 32 GeV 4 


CM 


1.3 x 10 16 GeV 


The flaton and .R-axion mass is estimated as 


m x ~ m a ~ 450 GeV 
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1.3 x 10 16 GeV' 
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5 x 10 " 8 
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(87) 


( 88 ) 


where the coupling constant |A^| is assumed to be small enough for the flaton mass to 
be smaller than sparticle mass, which keeps it from decaying into sparticle pairs. From 
Eq. (18), the VEV of the flaton at the true minimum is given by 

!/3 / I \ 1 \ -V3 

M ~ 8.6 x 10 13 GeV 1 ' 1 1 


CM 


1.3 x 10 16 GeV 2 
Then, the /i-terrn outside the Q-ball is given by 

I 


x\ 


5 x lO " 8 


Ml* out I) ~ 300 GeV 


10 


-7 


M 


8.6 x 10 13 GeV 


(89) 


(90) 


12 Since ( Vq/3M ^)F 2 > to^, the moduli field may start to oscillate after the potential energy of the 
flaton dominate the energy of the universe, which implies that the thermal inflation cannot work. Here, 
we assume that the moduli field starts to oscillate before the thermal inflation begins. 
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where the coupling constant X tl | is also assumed to be as small as 10“' to obtain the 
//-term of the electroweak scale. 

After the thermal inflation, the VEV of the flaton violates the condition for the exis¬ 
tence of the Q-ball and Q-balls decay. The released lepton number is converted into baryon 
number through the sphaleron process. The provided baryon asymmetry is estimated as 


Y b = 


~n L 


23 P V ,BB 


Pri, BB 


77 osc 


~ 7.7 x 10“ 10 


-RH 


m 3/ 2 

5 x 10 7 GeV ) V 200 keV 


-1 


x 




1.7 x 10 32 GeV 


rpX 

1 RH 


10 MeV 7 V 20 TeV 


T r 


M 


pi 


(91) 


One can hnd that the observed baryon asymmetry, Yb — 8 x 10“ n , could be explained if 
the AD held begins to oscillate near the Planck scale. Q-balls collapse irrespective of the 
oscillation amplitude of the AD held, which is contrary to the case of unstable Q-balls in 
Sec. 4.2. 


Next, we estimate the temperature at the end of the thermal inflation, T c , and the 
reheating temperature of the baton decay, T^ H . The thermal inflation ends when T c ~ 
c 1 / 2 (|X|) because thermal particles which couple with the baton become massive. From 
Eq. (21), the temperature at the end of the thermal inhation, T c , is estimated as 


T c ~ 17 TeV 


a 


1.3 x 10 16 GeV 2 


/ m 3/ 2 y / 3 
\ 200 keV / ’ 


(92) 


where we assume c ~ 0(1). Note that the released lepton number is successfully converted 
to baryon number through the sphaleron process since the thermal inhation ends before 
the electroweak symmetry breaking. As for the reheating temperature, T^ H is estimated 
from the decay rate of the baton. Since m x > 2m^, the baton mainly decays into two 
Higgs bosons at the tree level and they decay into the Standard Model particles. From 
Eq. (11), The reheating temperature, T^ H , is expressed as 



where we use g* = 10.8. Hence, the reheating temperature can be higher than 0( 10) MeV, 
which can avoid spoiling the success of the BBN. 
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We turn to estimate the density parameter of the moduli. The density parameter of 
the Big-Bang moduli after the entropy dilution is given by 


n^BB^ 2 — 6.1 X 10 5 


-EH 


5 x 10 7 GeV 


Vr 


0 


X 


7^X 

_ RH 

10 MeV ) \ 20 TeV 


1.7 x 10 32 GeV 4 
2 


-1 


Vo 


M, 


pi 


(94) 


It is found that the Big-Bang moduli is sufficiently diluted. On the other hand, the density 
parameter of the thermal inflation moduli is given by 


a 


■ 77 , TH 


~ 8.1 x 10 7 4 


On 


1.7 x 10 32 GeV 4 


m 3/2 \ ~ 2 
200 keW 


rpX 

_ RH 

10 MeV 



(95) 


One can find that the density of the thermal inflation moduli is much larger than the 
critical density. Therefore, in order for this scenario to work, the separation between the 
local minimum determined by the Hubble induced term and the true minimum should 
be of the order of 10 _4 M p i. 13 Although this may result from 0.01% fine-tuning of the 
moduli potential, we have no motivation for the thermal inflation in that case since the 
moduli problem can be solved when Trh — lOMeV and r] 0 /Mp\ ~ C>(1CD 3 ) (see Eq. (2)). 
Therefore, it is found that this scenario cannot work due to the secondary produced 
moduli. 


6 Summary and Discussions 

In gauge-mediated SUSY breaking models, the cosmic density of the coherently oscillating 
moduli exceeds the dark matter density unless it is diluted by some entropy production 
mechanism after the primordial inflation. The thermal inflation can successfully dilute the 
moduli abundance. However, there exists another fatal problem. The thermal inflation 
also dilutes baryon number which may be produced beforehand. Hence, some viable 
baryogenesis is needed for a cosmologically consistent scenario. A promising candidate 
for a mechanism to produce huge baryon asymmetry is the Affleck-Dine mechanism. In 
this paper, we have studied if the thermal inflation is consistent with the Affleck-Dine 
baryogenesis in gauge-mediated SUSY breaking models. 

13 If the moduli field has no couplings with photons, the constraint from the X-ray background spectra 
is irrelevant. In this case, the moduli field with mass of 0 ( 1 ) MeV, which cannot decay into electrons, 
could be the dark matter and the fine-tuning can be relaxed. However, it is still difficult for this scenario 
to work unless fine-tuning is allowed. 
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Q-ball formation is inevitable for the logarithmic potential induced from the gauge- 
mediated SUSY breaking effects and becomes a hinderance to the baryon number pro¬ 
duction. In order to produce huge baryon number, the amplitude of the AD held should 
be large, which renders the formed Q-balls stable. It is difficult to extract baryon charge 
from stable Q-balls and the evaporated baryon charge cannot explain the observed baryon 
asymmetry. On the other hand, unstable Q-balls can release all baryon charge. In this 
case, however, the amplitude of the AD held is restricted from above in order to prohibit 
Q-ball decay during and after the BBN epoch. Thus, sufficient baryon number cannot 
be produced. As a result, the Affleck-Dine mechanism is incompatible with the thermal 
inflation to solve the moduli problem in both cases. 

Among the hat directions, the LH U hat direction is special since it has the super- 
symmetric /i-term which violates the condition for the existence of the Q-ball. Since the 
baton can have a coupling with the the Higgs supermultiplets, the baton VEV naturally 
generates the fi- term. We showed that the generated /i-term at the end of the thermal 
inhation destroys formed Q-balls and that the released lepton asymmetry could explain 
the baryon asymmetry. In this case, however, the energy of the thermal inhation Vo is 
required to be large in order to violate the condition for the existence of the Q-ball, which 
leads to overproduction of the thermal inhation moduli. Thus, the moduli problem arises 
again because of the secondary produced moduli relic and serious hne-tuning is needed 
for this scenario to work. We conclude that it is difficult for the thermal inhation to be 
consistent with the Affleck-Dine mechanism. 

One might think that the clectroweak baryogenesis can explain the asymmetry without 
dilution. In that case, the reheating temperature after the thermal inhation, T^ H , should 
be higher than the electroweak scale. However, it seems to be difficult to dilute the moduli 
density for such high reheating temperature in the context of the thermal inhation. 

Finally, we make some comments on gravity-mediated SUSY breaking models. In 
these models, the gravitino mass is comparable to sparticle mass of the order of msusY — 
0(1) TeV. The usual Affleck-Dine mechanism cannot provide sufficient baryon number 
to survive after the dilution because of the early oscillation of the AD held by the soft 
SUSY breaking mass term. Instead, the modified Affleck-Dine mechanism, which has 
been proposed by [15], could explain the observed baryon asymmetry. In this model, the 
LH U hat direction provides lepton number after the thermal inhation. The dynamics of 
the LH U hat direction is so complicated that numerical simulations are necessary. Some 
works have revealed that this modified mechanism can work in the gravity-mediated SUSY 
breaking models [44, 45], but it is unclear if it can also work well in gauge-mediated SUSY 
breaking models. 
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